arXiv:1504.02034vl [physics.space-ph] 8 Apr 2015 


Monte Carlo model for electron degradation in methane gas 


Anil Bhardwaj*, Vrinda Mukundan 


Space Physics Laboratory, Vikram Sarabhai Space Centre, Trivandrum 695022, India 


Abstract 

We present a Monte Carlo model for degradation of 1-10,000 eV electrons in an atmosphere of methane. The electron 
impact cross sections for CH4 are compiled and analytical representations of these cross sections are used as input 
to the model.“Yield spectra”, which provides information about the number of inelastic events that have taken place 
in each energy bin, is used to calculate the yield (or population) of various inelastic processes. The numerical yield 
spectra, obtained from the Monte Carlo simulations, is represented analytically, thus generating the Analytical Yield 
Spectra (AYS). AYS is employed to obtain the mean energy per ion pair and efficiencies of various inelastic processes. 
Mean energy per ion pair for neutral CH4 is found to be 26 (27.8) eV at 10 (0.1) keV. Efficiency calculation showed that 
ionization is the dominant process at energies >50 eV, for which more than 50% of the incident electron energy is used. 
Above 25 eV, dissociation has an efficiency of ~27%. Below 10 eV, vibrational excitation dominates. Contribution 
of emission is around 1.2% at 10 keV. Efficiency of attachment process is ~0.1% at 8 eV and efficiency falls down 
to negligibly small values at energies greater than 15 eV. The efficiencies can be used to calculate volume production 
rate in planetary atmospheres by folding with electron production rate and integrating over energy. 
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1. Introduction 


Electron collision with molecule can result in various processes, like ionization, dissociation, and excitation of the 
target molecule, which can produce new species that can be more reactive as well as physically and chemically 
different from their parent molecule. Secondary electron released during ionization can also initiate further reactions. 
By-products of all these processes can initiate further reactions which are of great interest in the field of radiation 
chemistry, environmental chemis try, p lanetary aeronomy processes, like aurora and dayglow, and also in astrophysical 
and biological systems (Mason, 2003; Campbell & Brunger, 2009). To understand such phenomena, a thorough 
knowledge of electron degradation when it collides with atoms or molecules is required. 

Methane is the simplest hydrocarbon present in the solar system ( Mueller-Wodarg et all 2008 ). It causes infrared 
absorption in the atmosphere of Jupiter and Saturn and is an important atmospheric constituent in the planets Uranus 
and Neptune ( Broadfoot eUd. . 19791). In Titan, photochemistry is governed by ionization and dissociation products 
of nitrogen and methane dLavvas et al. . 2011 ). Collision of solar photons or photoelectrons with methane molecules 
causes the neutral dissociation or ionization of the molecule which in turn leads to the generation of simple hydro¬ 
carbon radicals and ions. The subsequent reactions caused by these radicals and ions, either with themselves or with 
methane and other background gases , cau se the production of higher order hydrocarbons, be it alkanes, alkenes or 
alkynes ( Banaszkiewicz etldl 200(1 Strobe! 1200 4) and leading to polymerization which may produce UV-dark haze 
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in auroral region of Jupiter ( Singhal et all 1992) and very heavy ionic species in Titan’s atmosphere ( Coates et all 
2007 ; Wahlund et al. , 20091) . Hydrogen cyanide, an important precursor f or the formation of am ino acids and proteins, 
is formed from those reactions for which methane acts as a precursor ( Romanzin et al. . 20051) . Fifth flyby of Titan 
by Cassini-Huygens mission found regions of low radar reflectivity which are interpreted as lakes, with methane as a 


major constituent (Cordiar et aj„ 2009). 


The aim of this study is to present a Monte Carlo model which describes the degradation of electrons with energy 
in the range 1 eV to 10 keV in a CH 4 atmosphere. Gan & Cravens (1992) used solution of Boltzmann equations 
for studying degradation of electrons in CH 4 and they calculated energy transfer rates for elastic and various in¬ 
elastic processes assuming a Maxwellian electron distribution. Monte Carlo method is a stochastic method, which 
has been w idely used for study i ng the problem of elec t ron energy degradation in gases relevan t for planetary at¬ 
mospheres llCicerone & Bowhill ( 1971 f Ashiharaj ( 1978b, Green et all ( J_977jk SinghaletalJj 1980b. Singhal_&_Green 


Singhal & Bha rdwaj (1 991 ) 


1981), 


( 2000 ), 


Bhardwaj & Micheal (1999b), Shematovich et al 


Bhardwaj & Singhal ( 


1993), Bhardwaj & Michael (1999a), Michael & Bhardwaj 


(2008), Bhardwaj & Jain (2009)]. In this method, history of 


collisions of particles is simulated, and conclusions are drawn from the statistics of those histories. Even though time 
consumin g, at so me levels it is found to be the most realistic simulations possible for studying electron energy depo¬ 
sition ( Solomon, 200lb . The energy loss process of electrons is actually discrete in nature and this nature is exactly 
captured in the Monte Carlo model. The method make use of probabilistic decision making techniques, and accuracy 
of the result largely depends on the number of simulations carried out. The study involves two steps: compilation of 
cross sections for all e-CH 4 collision processes and development of an energy apportionment method to determine 
how electron energy is distributed in various loss channels. 


2. Cross Sections 


2.1. Total and Differential Elastic 


Total el astic scattering c ross section for m ethane have been measured bv lBoesten & Tanakal (1 199 lb . iBundsch u et al 


(1997), Igaet al 


Measurements of 


19991) . and Kanik et al. (1993). All these measurements are in good agreement with each other. 


Boesten & Tanak a (1991) in the en ergy range 1 . 5-100 eV was fitted using analytical formula by 


Shirai et al. (2002). At energies above 100 eV, data of Kanik et al. (1993) has been used for fitting. This analytically 


fitted form of elastic cross section is used in the current study and is shown in Figure Q] 

The direction in which the electron is scattered after collision with a CH 4 molecule is determined using differential 
elastic cross sections (DCS). DCS for e-CH 4 collision has been measured by many workers. Values of DCS used in the 
presen t work are given in TableQ] In the low energy range of 3 to 15 eV, DCS measure ments of Manstone & Newell 
(11992b are used. However, DCS value at 5 and 10 eV are taken from lCho et alJ (12008b . Cross secti ons for energies 
between 20 to 100 eV also are taken from Cho et al. ( 2008 ). From 200 to 500 eV, measurements of Iga et al. (1999) 


and at 700 eV measurements of Sakae et al 


( 1989b are used. Since DCS measurements are not available for CH 4 for 


energies greater than 700 eV, linearly extrapolated values of cross sections are used. 
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Figure 1: Cross sections for elastic and inelastic processes for e-CFLt collisions. vi ,3 and V 2,4 are the cross sections 
for stretching and bending vibration modes, respectively. 
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Table 1: Elastic differential cross section for CH 4 in units of cm 2 . Value inside the bracket indicates a linearly extrapolated value. Notation IE-18 implies 1 
x 1 (T 18 


Angles 

n 

Energy 

(eV) 

0 

10 

20 

30 

40 

50 

60 

70 

80 

3.2 

(2.20E-17) 

(2.80E-17) 

(3.40E-17) 

4.00E-17 

4.60E-17 

5.40E-17 

9.40E-17 

8.10E-17 

8.50E-17 

4.2 

(5.70E-17) 

(6.40E-17) 

(7.10E-17) 

7.80E-17 

8.50E-17 

1.1 IE-16 

1.19E-16 

1.14E-16 

1.01E-16 

5 

(8.23E-16) 

(6.09E-16) 

3.91E-16 

1.75E-16 

1.32E-16 

1.28E-16 

1.39E-16 

1.48E-16 

1.39E-16 

6 

(3.19E-16) 

(2.73E-16) 

(2.27E-16) 

1.81E-16 

1.35E-16 

1.42E-16 

1.30E-16 

1.42E-16 

1.38E-16 

7.9 

(7.32E-16) 

(6.04E-16) 

(4.76E-16) 

3.48E-16 

2.20E-16 

1.80E-16 

1.43E-16 

1.39E-16 

1.32E-16 

10 

9.65E-16 

(8.03E-16) 

6.40E-16 

4.88E-16 

3.51E-16 

2.12E-16 

1.77E-16 

1.14E-16 

0.86E-16 

15.4 

(7.65E-16) 

(6.29E-16) 

(4.93E-16) 

3.57E-16 

2.21E-16 

1.75E-16 

1.08E-16 

6.50E-17 

5.00E-17 

20 

(1.20E-15) 

9.45E-16 

6.91E-16 

4.55E-16 

2.69E-16 

1.34E-16 

0.83E-16 

0.53E-16 

0.36E-16 

30 

(0.39E-15) 

10.07E-16 

6.25E-16 

3.36E-16 

1.60E-16 

0.75E-16 

0.40E-16 

0.27E-16 

0.18E-16 

50 

(1.45E-15) 

9.22E-16 

3.98E-16 

1.80E-16 

0.70E-16 

0.29E-16 

0.18E-16 

0.10E-16 

0.07E-16 

100 

(1.41E-15) 

8.01E-16 

1.95E-16 

0.44E-16 

0.19E-16 

0.09E-16 

0.04E-16 

0.03E-16 

0.02E-16 

200 

(1.14E-15) 

5.51E-16 

1.03E-16 

3.34E-17 

1.51E-17 

8.20E-18 

5.20E-18 

3.80E-18 

3.20E-18 

300 

(9.29E-16) 

4.26E-16 

7.26E-17 

2.56E-17 

1.12E-17 

6.20E-18 

3.80E-18 

2.60E-18 

1.80E-18 

400 

(6.76E-16) 

3.02E-16 

5.14E-17 

1.83E-17 

8.20E-18 

4.00E-18 

2.60E-18 

1.50E-18 

1.10E-18 

500 

(6.46E-16) 

2.80E-16 

4.58E-17 

1.63E-17 

6.90E-18 

3.20E-18 

1.80E-18 

1.20E-18 

9.00E-19 

700 

(1.54E-15) 

2.19E-16 

3.52E-17 

1.26E-17 

4.41E-18 

1.98E-18 

1.02E-18 

6.47E-19 

4.45E-19 

800 

(2.38E-15) 

(1.94E-16) 

(3.09E-17) 

(1.12E-17) 

(3.53E-18) 

(1.56E-18) 

(7.69E-19) 

(4.75E-19) 

(3.13E-19) 

900 

(3.67E-15) 

(1.71E-16) 

(2.71E-17) 

(9.74E-18) 

(2.82E-18) 

(1.22E-18) 

(5.78E-19) 

(3.49E-19) 

(2.20E-19) 

1000 

(5.67E-15) 

(1.5 IE-16) 

(2.37E-17) 

(8.56E-18) 

(2.25E-18) 

(9.64E-19) 

(4.35E-19) 

(2.56E-19) 

(1.55E-19) 

2000 

(2.56E-14) 

(6.66E-17) 

(9.8 IE-18) 

(3.66E-18) 

(5.09E-19) 

(2.04E-19) 

(6.70E-20) 

(3.33E-20) 

(1.54E-20) 

3000 

(4.56E-14) 

(4.12E-17) 

(5.85E-18) 

(2.22E-18) 

(2.13E-19) 

(8.27E-20) 

(2.24E-20) 

( 1 . 01 E- 20 ) 

(4.02E-21) 

4000 

(6.56E-14) 

(2.93E-17) 

(4.06E-18) 

(1.56E-18) 

(1.15E-19) 

(4.34E-20) 

(1.03E-20) 

(4.33E-21) 

(1.54E-21) 

5000 

(8.56E-14) 

(2.25E-17) 

(3.05E-18) 

(1.19E-18) 

(7.14E-20) 

(2.64E-20) 

(5.66E-21) 

(2.25E-21) 

(7.36E-22) 

6000 

(1.05E-13) 

(1.82E-17) 

(2.42E-18) 

(9.52E-19) 

(4.83E-20) 

(1.75E-20) 

(3.46E-21) 

(1.31E-21) 

(4.01E-22) 

7000 

(1.25E-13) 

(1.51E-17) 

(1.99E-18) 

(7.88E-19) 

(3.47E-20) 

(1.24E-20) 

(2.28E-21) 

(8.36E-22) 

(2.40E-22) 

8000 

(1.45E-13) 

(1.29E-17) 

(1.68E-18) 

(6.69E-19) 

(2.61E-20) 

(9.23E-21) 

(1.59E-21) 

(5.64E-22) 

(1.54E-22) 

9000 

(1.65E-13) 

(1.12E-17) 

(1.44E-18) 

(5.79E-19) 

(2.02E-20) 

(7.09E-21) 

(1.15E-21) 

(3.99E-22) 

(1.04E-22) 

10000 

1.85E-13 

(9.96E-18) 

(1.26E-18) 

(5.09E-19) 

(1.61E-20) 

(5.60E-21) 

(8.72E-22) 

(2.93E-22) 

(7.35E-23) 



Angles 

90 

100 

110 

120 

130 

140 

150 

160 

170 

180 

gree) 











Energy 











(eV) 











3.2 

8.30E-17 

6.20E-17 

3.40E-17 

2.20E-17 

2.00E-17 

2.70E-17 

(3.40E-17) 

(4.10E-17) 

(4.80E-17) 

(5.50E-17) 

4.2 

1.01E-16 

8.60E-17 

7.30E-17 

4.20E-17 

2.50E-17 

5.30E-17 

(8.10E-17) 

(1.09E-16) 

(1.37E-16) 

(1.65E-16) 

5 

1.19E-16 

0.84E-16 

0.50E-16 

0.24E-16 

0.18E-16 

0.26E-16 

0.36E-16 

0.47E-16 

0.60E-16 

0.70E-16 

6 

1.18E-16 

8.10E-17 

4.10E-17 

2.20E-17 

4.00E-17 

8.00E-17 

(1.20E-16) 

(1.60E-16) 

(2.00E-16) 

(2.40E-16) 

7.9 

9.80E-17 

6.20E-17 

2.70E-17 

1.70E-17 

4.60E-17 

9.40E-17 

(1.42E-16) 

(1.90E-16) 

(2.38E-16) 

(2.86E-16) 

10 

0.55E-16 

0.38E-16 

0.24E-16 

0.21E-16 

0.31E-16 

0.58E-16 

0.83E-16 

1.20E-16 

1.60E-16 

1.78E-16 

15.4 

4.00E-17 

2.90E-17 

2.70E-17 

3.20E-17 

4.30E-17 

5.70E-17 

(7.10E-17) 

(8.50E-17) 

(9.90E-17) 

(1.13E-16) 

20 

0.26E-16 

0.20E-16 

0.19E-16 

0.25E-16 

0.32E-16 

0.40E-16 

0.49E-16 

0.56E-16 

0.62E-16 

0.69E-16 

30 

0.12E-16 

0.12E-16 

0.13E-16 

0.17E-16 

0.21E-16 

0.27E-16 

0.30E-16 

0.32E-16 

0.33E-16 

0.34E-16 

50 

0.05E-16 

0.04E-16 

0.06E-16 

0.08E-16 

0.10E-16 

0.14E-16 

0.16E-16 

0.18E-16 

0.19E-16 

0.20E-16 

100 

0.02E-16 

0.02E-16 

0.03E-16 

0.03E-16 

0.04E-16 

0.04E-16 

0.05E-16 

0.05E-16 

0.05E-16 

0.05E-16 

200 

2.80E-18 

2.40E-18 

2.20E-18 

2.00E-18 

(2.15E-18) 

(2.25E-18) 

(2.35E-18) 

(2.45E-18) 

(2.55E-18) 

(2.65E-18) 

300 

1.40E-18 

1.20E-18 

1.10E-18 

1.10E-18 

(1.00E-18) 

(1.00E-18) 

(1.00E-18) 

(1.00E-18) 

(1.00E-18) 

(1.00E-18) 

400 

9.00E-19 

8.00E-19 

6.00E-19 

6.00E-19 

(6.00E-19) 

(6.00E-19) 

(6.00E-19) 

(6.00E-19) 

(6.00E-19) 

(6.00E-19) 

500 

7.00E-19 

6.00E-19 

5.00E-19 

5.00E-19 

4.00E-19 

4.00E-19 

4.00E-19 

4.00E-19 

4.00E-19 

4.00E-19 

700 

3.07E-19 

2.46E-19 

1.96E-19 

1.68E-19 

1.59E-19 

1.51E-19 

1.42E-19 

1.34E-19 

1.25E-19 

1.17E-19 

800 

(2.03E-19) 

(1.57E-19) 

(1.23E-19) 

(9.74E-20) 

(1.01E-19) 

(9.28E-20) 

(8.46E-20) 

(7.75E-20) 

(6.99E-20) 

(6.33E-20) 

900 

(1.34E-19) 

(1.01E-19) 

(7.68E-20) 

(5.64E-20) 

(6.32E-20) 

(5.70E-20) 

(5.04E-20) 

(4.49E-20) 

(3.91E-20) 

(3.42E-20) 

1000 

(8.92E-20) 

(6.46E-20) 

(4.81E-20) 

(3.27E-20) 

(3.98E-20) 

(3.50E-20) 

(3.00E-20) 

(2.59E-20) 

(2.18E-20) 

(1.85E-20) 

2000 

(6.13E-21) 

(3.41E-21) 

(2.21E-21) 

(9.06E-22) 

(1.89E-21) 

(1.41E-21) 

(9.88E-22) 

(6.93E-22) 

(4.66E-22) 

(3.24E-22) 

3000 

(1.28E-21) 

(6.1 IE-22) 

(3.65E-22) 

(1.11E-22) 

(3.20E-22) 

(2.16E-22) 

(1.34E-22) 

(8.35E-23) 

(4.93E-23) 

(3.05E-23) 

4000 

(4.21E-22) 

(1.80E-22) 

(1.01E-22) 

(2.51E-23) 

(9.06E-23) 

(5.71E-23) 

(3.25E-23) 

(1.85E-23) 

(1.00E-23) 

(5.70E-24) 

5000 

(1.78E-22) 

(7.00E-23) 

(3.78E-23) 

(7.91E-24) 

(3.40E-23) 

(2.03E-23) 

(1.08E-23) 

(5.79E-24) 

(2.90E-24) 

(1.55E-24) 

6000 

(8.80E-23) 

(3.23E-23) 

(1.68E-23) 

(3.08E-24) 

(1.52E-23) 

(8.75E-24) 

(4.42E-24) 

(2.23E-24) 

(1.05E-24) 

(5.35E-25) 

7000 

(4.85E-23) 

(1.68E-23) 

(8.49E-24) 

(1.38E-24) 

(7.77E-24) 

(4.28E-24) 

(2.06E-24) 

(1.00E-24) 

(4.49E-25) 

(2.18E-25) 

8000 

(2.89E-23) 

(9.54E-24) 

(4.69E-24) 

(6.95E-25) 

(4.32E-24) 

(2.31E-24) 

(1.07E-24) 

(4.98E-25) 

(2.14E-25) 

(1.00E-25) 

9000 

(1.83E-23) 

(5.78E-24) 

(2.78E-24) 

(3.78E-25) 

(2.58E-24) 

(1.33E-24) 

(6.00E-25) 

(2.69E-25) 

(1.11E-25) 

(5.03E-26) 

10000 

(1.22E-23) 

(3.70E-24) 

(1.74E-24) 

(2.19E-25) 

(1.62E-24) 

(8.22E-25) 

(3.57E-25) 

(1.55E-25) 

(6.21E-26) 

(2.72E-26) 




Figure 2: Electron impact cross sections of CH 4 for various ionization and dissociation channels. 


2.2. Attachment 


Dissociative electron attachment process of CH 4 results in the production of H and CfL ions. Cross section for this 
process was measured by Sharp & Dowell ( 19671) and Rawat et al. ( 2007 ). Former cross sections were analytically 
fitted by Shirai et al. (2Q02), which are used in the present work and is shown in FigureQ] 


2.5. Vibrational Excitation 


Methane molecule is found to have four normal modes of vibration: vi with threshold energy 0.362 eV, V 2 with 0.190 
eV, V 3 with 0.374 eV and V 4 with 0.162 eV. But it is difficult to resolve these modes experimentally as they have very 
close transition energies. Hence, experimental data are available for the combined cross section for symmetric vi 
and antisymmetric V3 stretching vibrations (vi,3), and symmetric V2 and antisymmetric V4 bending vibrations (>2,4) of 
CH 4 . 


Shvnl (1991) measured vibrational excitation cross sections for methane at 5.0, 7.5, 10.0 and 15 eV. Tawara (1992) 
mea sured cross sections in the energy range 0.16-100 eV. DCS values in the 0.6 to 5.4 eV range, measured by 


Bundschu et al. (1997), were integrated to obtain the integral cross section. Vibrational excitation cross sections 


used in the current study are taken from Davies et al. (1988) in which cross section values are given for a larger en¬ 
ergy range; 0.450 to 100 eV for vi ,3 mode, and 0.162 to 150 eV for V 2,4 mode. The values of Davies et al. (1988) 
and Tawara ( 1992h agree well at energies 10-100 eV. However, at energies less than 10 eV there is a difference be- 
tween two cross se ctions to a maximum of 5 0% at few energ i es. G ood agreement is f ound when c ross sections of 
Davies et al. ( 1988 ) are compared wi t h that of Bundschu et al. ( 1997 ). Measurements of Shvn ( 199 lh are found to be 


lower than the values of Davies et al. ( 19881) by ~50%. Vibrational cross section used in our study is shown in Figure 

□ 


2.4. Ionization 


Ionization and dissociative ionization of CH 4 results in the production of ions CH 4 + , CH 3 + , CH 2 % CH + , CV, H 2 _and 
H + . Cross sections for these ionization proce sses have been measured by many authors, e .g., iTian & Vidall 0997), 


Chatham et al. ( 1984 ). Adamczvk et al. d 19661) and Straub et~akl ( 1997 ). Straub et ah ( 1997 ) have measured the 


cross 


sections in the energy range 15-1000 eV and it is found to be the most reliable among various available measure¬ 
ments (Liu & Shemanskv , 20061) . These measurements were later revised due to instrumental recalibration and was 
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Table 2: Parameters for CH4 + , CH 3 + and CH2 + ionization cross section (equation (03 ) 


Process 

Ert(eV) 

ai 

a2 

a 3 

a 4 

as 

ae 

ch 4 + 

12.99 

4.40 

1.627 

7.720E-3 

-4.50E-2 

3.10E-2 

0.93 

ch 3 + 

14.24 

2.18 

1.435 

1.13E-2 

7.4E-2 

4.91E-2 

1.01 

ch 2 + 

15.20 

0.121 

1.868 

3.44E-2 

3.00E-1 

5.20E-2 

0.91 


publis hed in Lindsay & Manganl (20031 (Here after referred to as Straub’s revised measurements). Liu & Shemanskv 
(2006) derived the oscillator strength and excitation functions for various ionization channels of CH4 and calculated 
the cross section values. These cross sections are in good agreement with Straub’s revised measurements. However 
cross sections reported by |Liu_& Shemanskv! ( 20061) exclude the contribution by pair production (e.g.(CH2 + , H + ), 
(C + ,H + )). Erwin & Kune (12008 ) using scaling law, developed analytical expressions for calculating cross sections 
for various ionization channels of methane, which are valid at all non-relativistic energies. These expressions allow 
calculation of the electron impact ionization cross sections in an easier, and more direct way than the functions derived 
by Liu & Shemanskv ( 20061) . Good agreement is found between these theoretical ionization cross sections when com¬ 
pared with Straub’s revised measurements and measurements of Tian & Vidal ( 1997t) . Maximum deviation (-20%) is 
found for H + production. 


For the present work, we have taken the cross sections for CH4 
For CH 2 


and CH 3 + production from Liu & Shemanskv (2006). 
production, the cross section for positive ion pa i r form ation (CH2 + , H + ), measured bv lLindsav et alJ (12001). 


have been added with the values of Liu & Shemanskv! ( 20061) to account for the c ontribution via doubly ionized 
channels. These cross sections are fitted using analytical equation ( Shirai et al. . 20021) : 


cr = <r 0 a x (EIE R T-l [1 + (£/a 3 ) fl2+fl4 + (E/a 5 y 2+a6 ] 


(1) 


where cr 0 — 1 x 10 16 cm 2 , E« is the Rydberg constant, and ai, a?, a 3 , a 4, as and a () are the fitting parameters whose 
values are given in Table[2] The analytical expression of Erwin & Kune ( 120081) is used to calculate the cross sections 
for channels CH + , C + and FL \ For H + production channel, we have used the analytical expression of 1 Shirai et al 


(2002), extending it to 10 keV to get cross section values at higher energies. Figure 0 shows the ionization and 
dissociative ionization cross sections used in the present study. 


2.5. Dissociation 


Dissociation of methane by electron impact results in the production of neutral radicals CH 3 , CH 3 and CH. Exper¬ 
imental cross sections for these processes are not available over a wide energy range , except for the CH 3 r adical 
production where the measurements are made up to 500 eV ( Motlagh & Moore , 1998 ). Erwin & Kune ( 200^ ) have 
given analytical expression for CH4 dissociation cross sections which is valid at all non-relativistic energies. But 
the analytical expression does not account for the production of CH 3 radical through disso ciati ve ionization channel 
CH 3 + H + . For CH 3 radical production, we have used the analytical representation of Motlagh & Moore : ( 1998 ) cross 
sections, as given by Shirai et all (2002) and ex tended it to 10 keV. Fo r CH2 and CH radical production, cross section 
are calculated using the analytical expression of Erwin & Kune ( 2008 ). Figure[2]shows the dissociation cross sections 
used in the study. 


2.6. Emission 


Electronically excited state of neutral CH4 leads to the dissociation of the molecule resulting in the production of ex¬ 
cited fragments ( Danko e t all 201 lh. Cross s ections for hydrogen lym an series and carbon lines for energies less than 
400 eV was measured by Pang et al~! (1 9871 ). Motohashi et al. (1 996 ) measured cross sections for the emission from 
various excited fragments : the hydrogen Lyman and Balmer series, CH band emission at 420-440 nm, line emissions 
from C at 165.7 nm and 156.1 nm and these cross sections are found to be in agreement with that of Pang et ak (1987). 
The uncertainties in these cross sections were estimated to be ±20% for H Lyman-a, ±12% for H Balmer-a, ±20% 
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Figure 3: Electron impact cross sections for H Lyman and H Balmer emissions. 


for CH band, and +50% for atomic carbon emission. Maximum energy of experimentally measured cross sections are 
1 keV, 6 keV, 5 keV and 1 keV for H Lyman-a-, H Balmer-ar, CH band, and C emissions, respectively. Analytic cross 
sections for these emission processes given by Shirai et al. (2002) are extended upto 10 keV and used in the current 
model, (cf. Figure [3} 


The total inelastic cross section is obtained by adding the cross sections of above mentioned inelastic processes. Cross 
sections for the emission process and H + production channel are not taken into account while calculating total inelastic 
cross section as they are already accounted in other channels. Figures[l] |2l [3] and[4]shows the cross sections for elastic 
and various inelastic processes of methane that have been used in the present work. 


3. Monte Carlo Model 

A model for local degradation of electrons in CH 4 gas is developed in the energy range 1 eV to 10 keV using the 
Monte Carlo technique. Electrons incident upon CH4 molecules deposit their energy into the gas through elastic or 
inelastic collisions. Inelastic collisions lead to further ionization, dissociation, excitation or attachment processes. 
Thus, there exists a range of possible channels through which an electron can degrade its energy on colliding with 
molecules. In the model, every incident electron is followed in a collision-by-collision manner until its energy goes 
below an assigned value. 

The Monte Carlo simulation starts after the initial energy of electron is assigned. The direction of the incoming 
particle (9, ft) is decided using random numbers Ri and R 2 as 


9 = cos _1 (l -2Ri), 

<p — 2nR 2 . 

The distance to next collision is calculated from 

S = — log( 1 - Ri)/ncTT , 


( 2 ) 

(3) 

(4) 
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Figure 4: Electron impact cross sections of CH 4 for CH band and various C 1 line emissions. 


where R 3 also is a random number and n is the number density of the target particles, which is taken as 10 10 cm 
crj is the total electron impact collision cross section, which is 


Ct - vel + & in 


(5) 


where cr e i and <x,„ are the total elastic and inelastic collision cross sections. Decision on the type of collision that occur 
is made by comparing the probabilities of elastic and inelastic collisions, P e / and P,„, with another random number 
generated, R 4 , where P,./ and P,„ are calculated as cr e // ir-/ and cr in / cr T . If P,./ > R 4 , an elastic collision has taken place. 

Energy loss that occurs during elastic collision, due to target recoil, is calculated as 

m 2 v 2 nrvVi cos <5 

A E= ---, (6) 

m + M m + M 


Vi = v 


m cosd [M 2 + m 2 (cos 6 - l )] 1 / 2 
in + M 111 + M 


Here 6 is the scattering angle in the laboratory frame, v and m are, respectively, the velocity and mass of the electron, 
and M is the mass of the target particle. Differential elastic cross sections (discussed in section 2.1) are used to obtain 
the scattering angle <5. Differential cross sections are fed numerically in the Monte Carlo model at 28 unequally spaced 
energy points (3.2, 4.2, 5, 6 , 7.9, 10, 15.4, 20, 30, 50, 100, 200, 300, 400, 500, 700, 800, 900, 1000 eV; and 1, 2, 3, 4, 
5, 6 , 7, 8 , 9, 10 keV) and at 19 scattering angles (0°, 10°, 20°, 30°, 40°, 50°, 60°, 70°, 80°, 90°, 100°, 110°, 120°, 130°, 
140°, 150°, 160°, 170°, and 180°). At intermediate energies and scattering angles, values are obtained through linear 
interpolation. The energy A E is subtracted from the energy of the primary particle. After the collision, the deflection 
angle relative to the direction (9, tp ) is obtained as 


cost? -cos 9 cos 9 - sin 9 sin 9 cos (p , 

cos (p — (cost? cos </> sin 0 sin0 - sin <f> sin 9 sin0 + sin t?cos 0cos 9 )/ sin 6 , (7) 

sin 0 = (cos t? cos 0 sin f? cos</> -cos^sint? sin</> + sin 9 sin <p cos 9 )/ sin 9 . 

Here 9 , (p are the scattering angles. 
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Table 3: Parameters for Secondary Electron Energy 


Process 

I (eV) 

Ts 

Tb 

T s 

Ta 

T b 

ch 4 + 

12.99 

30 

10 

-0.989 

10 

44.52 

ch 3 + 

14.24 

23 

28 

0.947 

100 

64.52 

ch 2 + 

15.2 

7.52 

-38.1 

2.6 

21.8 

44 

CH + 

22.6 

11 

-0.5 

-0.8473 

0 

0 

C + 

27 

11 

-0.5 

-1.8473 

200 

200 

H 2 + 

23.53 

2 

0.0 

-5.8473 

0 

0 


If P e i < R 4 and 
inelastic event, 
electron is then 


P,„ > R 4 , an inelastic collision has taken place. In this case decision has to be 
If it is an ionization event, then a secondary electron will be generated. The 
calculated with a random number R using the relation (Bhardwai & Jain, 2009) 


made on the type of 
energy of secondary 


T = I" ^ [tanf/e/T, +(R- l)/f 2 )] + T s 

e v + r B 


Ta 

E v + Tg 


( 8 ) 


where 


K\ = tan 


(E v - 1) 


-Ts + 


Ta 


(E v + T b ) 


r,y Ey 
(Ey + Tfi) 


Ki — tan 


Ts - 


T a 


(Ey + T b ) 


T S Ey 
(Ey + Tfi) 


Here E v is the incident electron energy, IT,, Ta, Ta , Tb, and T s are the fitting parameters which are given in Table 
[3 and I is the ionization threshold. Secondary electron energy, calculated using equation ([ 8 ]). is compared with the 
lower cut off energy, which is taken as 1 eV in the current model. If the energy of secondary electron is found to be 
greater than the cut off energy, then they are also followed in the same manner as that of a primary electron. Similarly 
the tertiary, quaternary, etc electrons are also followed in the Monte Carlo simulation. 


Every single electron is followed in similar manner and number of collision events occurred are recorded in appro¬ 
priate energy bins. The energy bin size is taken as 1 eV for the entire energy range. After each event, the amount of 
energy loss due to the event is subtracted from the electron energy. After subtraction, if the electron energy is higher 
than cut off energy, it is again followed in the simulation. Electron degradation continues until all energies are below 
1 eV. The sample size is taken as 10 6 for each simulation. 


4. Yield Spectrum 


A two dimensional yield spectrum, U(E,E 0 ), which is a function of spectral energy E and incident electron energy E c , 
is obtained as the output of the Monte Carlo simulation. Yield spectrum gives information about the number of energ 
loss events that have occurred in an energy bin and is defined as ( Bhard wai & Michael, 1999a; Bhardwai & Jai 


a 


200 91 


U(E,E 0 ) = 


ME) 
A E ’ 


(9) 


where N(E) is the number of inelastic collision events for which the spectral energy of the electron is between E and 
E + A E, with A E being the energy bin width, which is 1 eV in the current model. Figure [5] shows yield spectrum at 
five incident energies. If vibrational cross section measurements of Tawaral ( 19921 ) are used in the energy range 1-10 
eV, there is no significant change in yield spectrum. The maximum deviation is about to be 4.5%. 


10 







































For practical applications, yield spectrum is usually represented in the form 


U(E, E 0 ) = U a (E, E 0 ) H(Eq -E- E m ) + 6(E 0 - E). 


(10) 


Here H is the Heavyside function, with E m being the minimum threshold of the processes considered, and 8{Eq - E) 
is the Dirac delta function which accounts for the collision at source energy Eo. In many atmospheric applications, it 


is convenient to represent the yield spectrum, U a (E, Eo), in an analytical form (Green et aj„ 1977), 


U a (E,E 0 ) = A l { s 0 +A 2 (£ 1 0 - t /e 3l2+r ) 


( 11 ) 


where f = Eo/lOOO and e = E/I (/ is the lowest ionization threshold which is equal to 12.99 eV), Aj = 0.024, A 2 = 
4.40, t — 0, r — -0.27, and 5 = -0.085 are the best fit parameters. 


The yield spectra at energy region very close to E 0 shows rapid oscillation which is caused by Lewis effect. This 
irregular nature occurs due to the fact that the process of energy degradation is not continuous, but discrete in nature. 
There are only certain energies near E 0 , which an electron can acquire. For a process with threshold E m , the electron 
will suffer a minimum energy loss of E 0 - E,„ and this will bring down the electron energy to specific discrete values. 
No energy value in between E 0 - E,„ can be acquired by the electron. To account for this effect, Heavyside function is 
included in the right hand side of equation (ITOt . 


The analytical yield spectrum (AYS), given by equation (fill) , well represent the numerical yield spectrum at energy 
values greater than the ionization threshold (13 eV). To improve AYS at lower energies, we have used the additional 
function introduced by Bhardwai & Jain (l2009if . 


U b (E,E 0 ) 


EoAoe- A ’ x /A 3 
(1 + e AiX ) 2 


( 12 ) 


Here x = (E-A^/A^, and Aq, A 3 M 4 , A 5 and A b are the fitting parameters. The values are Ao = 0.9, A 3 = 8 . 5 ,A 4 = 7.0, 
A 5 = 0.001 and A () = 4.5. We have introduced the parameters A 5 and A () to get a better fit in the lower energy (<10 
eV) region. The final AYS is the sum of equations (ITTb and ( 1 1 21 . The numerical yield spectrum obtained from the 
model as well as analytical yield spectrum are shown in Figure 0 

The analytical property of AYS is very useful in determining various property of the gas, like mean energy per ion 
pair and efficiency. The population or yield of any state j, which is the number of inelastic events of type j caused by 
an electron while degrading it’s energy from E 0 to cut off, can be calculated using AYS as 

rEo 

Jj(Eo) = U(E,E Q )Pj(E)dE. (13) 

Jw, h 

Here W t h is the threshold for the j'* event and P j(E) is it’s probability at the energy E, which is calculated as P/E) = 
crj(E)/o-j„(E); cr in (E) being the total inelastic collision cross section at energy E. Yield of any process calculated using 
equation ([13} can be used to obtain mean energy per ion pair and efficiency. 


5. Mean Energy Per Ion Pair 

Mean energy per ion pair is defined as the average energy spent by an electron to produce an electron-ion pair after 
its energy is completely dissipated. Its reciprocal gives the efficiency with which a particle can ionize the gas, and is 
a typical feature of the target species considered. It is calculated as 


Ej(Eo) = E 0 /Jj(E 0 ), (14) 

where Jj(Eo) is the population of the f h process at the incident electron energy E 0 . At high incident electron energies, 
li approaches a constant value. Figure [ 6 ] shows /1 value calculated for neutral CH 4 and for the various ionization 
channels of methane. At ionization threshold, /1 shows a very high value. As incident electron energy increases. 
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Figure 5: Yield spectra for different incident energies. Solid curve shows numerical yield spectra obtained using the 
model. Analytical Yield Spectrum (AYS), calculated using equation (fTTh . is represented by dashed curves. Dotted 
curve shows improved AYS, obtained by summing equations (fTTb and (IT2l> . Yield Spectrum for 10000, 5000, 1000, 
100 and 50 eV are shown after multiplying with scaling factors 9000, 1000, 100 and 30 and 1, respectively 



Figure 6 : Mean energy per ion pair for ions CH 4 A CH 3 + , CH 2 + , CH + , C + , Hi + , H + , and neutral CH 4 (shown as total). 
Symbols shows the values calculated using numerical yield spectrum for CH 4 + and neutral CH 4 
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Energy (eV) 

Figure 7: Energy Distribution of secondary electrons for incident energies 100 eV, 1 keV and 10 keV. Y axis shows the 
number of secondary electrons produced per incident primary electron. Dotted and dashed curve shows distribution 
of tertiary and quaternary electrons, respectively, for an incident energy of 1000 eV. The inset shows the energy 
distribution of secondary, teritary and quaternary electrons for an incident energy of 10 keV by zooming in the lower 
energy range. 


population of ionization process increases as a result of which /j falls off rapidly. From about 100 eV onwards, the 
curve falls off very slowly and attains a constant value at high incident energies. The value of ft for CH 4 + , CH 3 + , 
CH 2 + , CH + , C + , H 2 + and H + ions are, respectively, 51.2 (54.9) eV, 68.3 (75.8) eV, 357.09 (369.9) eV, 972.7 (855.7) 
eV, 3.2 (2.7) keV, 4.08 (3.1) keV and 413.06 (358.16) eV, at an incident energy of 10000 (100) eV. The mean energy 
per ion pair for neutral CH 4 is 26 eV at 10 keV and 27.8 eV at 100 eV. Experimentally determined value for mean 

i ICRUl ( 1 9931) . is 27.3±0.3eV for incident energies > 10 keV. Fox et al. ( 20081) reported 


energy per ion pair, as given in 


a value of 31 eV, while Wedlund et al. 


( 201 lh computed a value of 28.0+1.2 eV at an incident electron energy of 2 


keV. Our value of mean energy per ion pair is consistent with those reported in previous studies. 


6. Secondary Electron production 

The secondary electrons that are produced during ionization events can have a maximum energy of (E-I)/2, where I is 
the ionization threshold. The energy of these electrons are calculated using equation ©. If the energy of the secondary 
electron is greater than that of the cut off energy, then it is also followed in the same manner as that of the primary. 
Similarly, tertiary, quaternary, etc electrons are also followed in the Monte Carlo simulation. The energy distribution 
of secondary electrons is shown in Figure |7]at few incident energies. Distribution of tertiary and quaternary electrons 
are also shown for incident energy of 10 keV. Figure [7] shows that, each incident electron of energy 10 keV, at some 
point of its energy degradation process, produces at least one secondary, or teritary or quaternary electrons whose 
energy is <32 eV, which is sufficient enough to cause further inelastic collisions. 


7. Efficiency 

During the degradation process, the electron energy is distributed among various inelastic processes. Efficiency of a 
process gives information on what fraction of the incident energy is used for a particular process after the electron has 
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Figure 8: Efficiencies of various ionization processes. Symbols represent the efficiencies that are calculated using 
numerical yield spectra and solid lines are efficiencies calculated using AYS. 


completely degraded its energy. The efficiency, rjj(Eo), of the f h process at the incident energy Eq can be obtained as 

rij(E 0 ) = Jj(E 0 ) (15) 

Eq 

where W t h is the threshold for the j' h process. The efficiency of various processes is calculated using numerical yield 
spectrum as well as AYS. 

Efficiencies of various ionization processes are shown in Figure [8] Because of its higher cross section, CH 4 + pro¬ 
duction channel has the highest efficiency throughout the energy range, with an efficiency of 25.3% (23.6%) for an 
incident electron energy of 10 keV (100 eV). The other ionization channels CH 3 + , CH 2 + , CH + , C + , H 2 + and H + have 
efficiencies of 20.8% (18.8%), 4.3% (4.1%), 2.5% (2.8%), 0.9% (1%), 0.5% (0.6%) and 4.4% (5%), respectively. At 
electron energies >100 eV, there is no significant variation in these efficiencies. But at lower energies, especially near 
the threshold region, ionization efficiencies fall off very rapidly. At 18 eV, the efficiencies for the production of CH 4 + , 
CH 3 + and CH 2 + are 18.2%, 9.1% and 0.2%, respectively. 

Figure[9]show efficiencies of various dissociation channels. Since the production of CH; radical has the highest cross 
section (cf. Figure [2}, it has the highest efficiency with a value of 20.8% (21%) at 10 keV (100 eV). Efficiencies 
of CFE and CH production are 3.9% (3.7%) and 2.5%(2.6%) at 10 keV (100 eV). The dissociation efficiency is 
almost constant at energies >100 eV. At electron energy of 30 eV, CH 3 , CH 3 and CH dissociation channels are having 
efficiencies 23.8%, 2.9% and 1.3%, respectively. 

Efficiencies of various emission processes are shown in Figures [TO] and [IT] Only a small fraction of incident electron 
energy goes to various emission channels with H Lyman-a emission having the highest efficiency of 0.43% (0.58%) 
at 10 keV (100 eV). For H Lyman-/1 and Lyman-y emissions, efficiencies are 0.11% (0.14%) and 0.05% (0.06%) at 10 
keV (100 eV). The CH band emission has an efficiency of 0.25% (0.28%) at 10 keV (100 eV). Among the various line 
emissions of atomic carbon, the 165.7 and 156.1 nm emission have almost the same efficiencies at all energies, with 
a value of 0.033% (0.04%) and 0.031% (0.03%) at 10 keV (100 eV), respectively. The carbon 193.1 nm emission has 
an efficiency of 0.02% at 10 keV. As there is a large uncertainty in the value of C-line emission cross sections (+50%), 
the calculated value of efficiencies would also be uncertain by similar amount. 

In Figure[12]an overall picture of efficiencies of various inelastic loss processes is presented. Efficiency values calcu¬ 
lated using both the numerical yield spectrum as well as the AYS are shown; a good match is observed between the 
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Figure 9: Efficiencies of various dissociation channels. Symbols represent the efficiencies that are calculated using 
numerical yield spectra and solid lines are efficiencies calculated using AYS. 



Figure 10: Efficiencies of various emission channels. Symbols represent the efficiencies that are calculated using 
numerical yield spectra and solid lines are efficiencies calculated using AYS. 
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Figure 11: Efficiencies of CH and various C emissions. Symbols represent the efficiencies that are calculated using 
numerical yield spectra and solid lines are efficiencies calculated using AYS. 


two efficiency values at energies greater than 10 eV. Efficiency calculated using the AYS for energy <10 eV would 
be quite approximate as the AYS is not able to represent well the numerical yield spectrum in this region. Hence, 
the efficiency of vibration process shown in Figure[T2]is calculated only using the numerical yield spectrum. Among 
the different loss processes ionization is found to be the dominant process above 30 eV. Above 100 eV, the ionization 
efficiency attains a constant value of 54%. The dissociation efficiency is constant at energies above 30 eV with a 
value of 27%. The emission efficiency is 1.2% (1.6%) at 10 keV (100 eV). All through the energy range, only a very 
small fraction of the incident electron energy is channeled into attachment process. The attachment efficiency peaks 
at 10 eV and has a value of 0.14%. 


8. Summary 


We have developed a Monte Carlo model for studying the degradation of 1-10,000 eV electrons in methane gas. 
Analytically fitted cross sections are used as input to the model. The numerical yield spectra (NYS) obtained as 
the output of the Monte Carlo model incl udes non-spatial in form ation about the degradat ion process. The NYS is 
analytically fitted using equations given by Green et al. (1977) and Bhardwai & Jain (20091). We have introduced two 
new parameters to better fit the low energy (< 10 eV) region of the NYS as described in equation[T2] thus obtaining 
the Analytical Yield Spectra (AYS). The AYS is used to calculate various parameters, like mean energy per ion pair 
and efficiency of various loss channels. The mean energy per ion pair for CH 4 has a value 26 (27.8) eV at 10 (0.1) 
keV. The energy distribution of secondary electrons for a few incident energies is presented in Figure |7| 


Efficiency of a loss channel gives information on the amount of incident electron energy going into that loss process. 
Efficiencies are calculated using the AYS as well as the NYS, and are found to be in good agreement for energy >10 
eV. At energies <10 eV, vibration is the dominant loss process with an efficiency 80% at 8 eV. In this energy region, 
electron attachment process has an efficiency of 0.14%, which falls down to very small value for energy >20 eV. 
From 25 eV onwards, dissociation process has an efficiency of 27%. At energies higher than 100 eV, ionization is the 
dominant loss process consuming more than 50% of the incident electron energy. In addition to the major inelastic 
processes, efficiencies are calculated for various emissions. The H Lyman-or emission has the highest efficiency among 
various emission channels. 

The results presented in the paper will be useful for modeling of aeronomical processes in the planetary atmospheres 
where methane is a significant constituent. Using AYS, photoelectron fluxes in the atmosphere can be calculated 
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Figure 12: Efficiencies of various loss channels. Symbols represent the efficiencies that are calculated using numer¬ 
ical yield spectra and solid lines are efficiencies calculated using AYS. Vibrational efficiencies are calculated using 
numerical yield spectrum only. 


which can be employed later on for calculating electron impact excitation or emission rates dJain & BhardwaiL 1201 1 


Bhardwai & Micheal , J_999bj). Energ y deposition rate can be calculated as a product of ionization rate and mean 


energy per ion pair ( Fox et all 2008 ). These efficiencies can be applied to planetary atmospheres for calculating 
volume production rates by multiplying with electron production rate and integrating over energy 
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